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AN  ANALOG  STUDY  OF  THE  RELATIVE  IMPORTANCE  OF 
VARIOUS  FACTORS  AFFECTING  ROLL  COUPLING 
By  Joseph  Weil  and.  Richard.  E.  Day 


SUMMARY 


An  analog  study  of  the  roll-coupling  problem  has  been  made  for  a 
representative  swept -wing  and  tailless  delta-wing  configuration.  The 
investigation,  made  primarily  for  subsonic  flight  conditions,  included 
the  determination  of  the  effects  of  wide  variations  in  many  of  the  per¬ 
tinent  aerodynamic  derivatives  on  the  motions  developed  in  rolling 
maneuvers.  The  influence  of  large  changes  in  principal  axis  inclination 
and  mass  distribution  is  also  considered. 

The  results  indicated  that  as  first  predicted  in  MCA  TN  1627  the 
relationship  between  the  longitudinal  and  directional  stability  is  of 
paramount  importance.  For  most  current  designs  an  optimum  condition 
exists  when  the  natural  frequencies  in  pitch  and  yaw  are  approximately 
equal.  Increases  in  pitch  damping  had  a  pronounced  favorable  effect  in 
reducing  the  amplitudes  of  the  motions  encountered  and  were,  in  general, 
considerably  more  effective  than  corresponding  increases  in  yaw  damping. 

Practical  redistribution  of  mass  produced  only  relatively  minor 
changes  in  the  overall  results. 

It  was  found  that  the  amplitude  of  the  motions  developed  for  a  given 
aileron  deflection  depends  to  a  large  extent  on  the  duration  of  the 
maneuver  (change  in  bank  angle).  limited  studies  indicated  that  90°  roll 
maneuvers  would  be  considerably  less  critical  than  j6o°  rolls.  The  angle 
of  attack  of  the  principal  axis  has  an  important  bearing  on  the  behavior, 
particularly  in  the  absence  of  other  disturbing  functions.  If  the  ini¬ 
tial  angle  of  attack  is  maintained  constant,  a  reduction  in  altitude  will 
delay  critical  conditions  to  a  higher  roll  rate  but  the  maximum  ampli¬ 
tudes  may  be  only  slightly  affected. 

Small  inadvertent  stabilizer  inputs  can  greatly  affect  the  magni¬ 
tude  of  the  motions  developed. 

It  is  difficult  to  generalize  on  the  effects  of  Mach  number  varia¬ 
tion  because  this  variable  affects  many  of  the  controlling  parameters. 
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Utilizing  simple  concepts  proved  useful  in  assessing  the  qualitative 
effects  of  many  of  the  aerodynamic  and  inertia  parameters  and  changes  in 
flight  condition.  The  calculated  lower  resonant  frequency  generally  cor¬ 
responded  to  the  average  roll  velocity  at  which  the  more  serious  motions 
could  he  expected. 

It  is  indicated  that  rational  design  procedure  can  avoid  the  problem 
of  serious  roll  coupling  at  supersonic  speeds  and  minimize  the  problem 
in  the  subsonic  speed  range  for  the  configurations  of  the  type  considered. 


INTRODUCTION 


Until  recently'the  motions  of  a  rigid  airplane  with  deflected  ailerons 
could  be  adequately  predicted  in  most  instances  from  the  classic  three- 
degree-of -lateral  freedom  solutions  (for  example,  ref.  l).  That  the 
longitudinal  modes  can  exert  a  powerful  effect  on  the  overall  motions  of 
the  rolling  aircraft  was  first  demonstrated  theoretically  in  19^8  (refs.  2 
and  3)i  however,  not  until  195^  were  strong  manifestations  of  coupled 
motions  in  rolling  maneuvers  experienced  on  a  full-scale  airplane  (ref.  4). 

Following  the  flight  experiences  described  in  reference  4,  a  five- 
degree-of -freedom  analog  study  was  initiated  by_  staff  members  of  the 
NACA  Langley  Laboratory  and  NACA  High-Speed  Flight  Station  to  determine 
whether  the  actual  flight  motions  could  be  calculated  and  also  to  deter¬ 
mine  the  effects  of  variations  in  certain  of  the  stability  parameters. 

Some  of  the  results  of  these  preliminary  studies  were  reported  briefly 
in  references  5  and  8. 

An  analog  investigation  of  several  generalized  airplane  configura¬ 
tions  encompassing  a  considerably  broader  scope  than  the  work  of  refer¬ 
ences  5  and  8  has  recently  been  completed  at  the  NACA  High-Speed  Flight 
Station.  In  this  study  wide  variations  in  many  of  the  pertinent  aerody¬ 
namic  derivatives  were  investigated  at  subsonic  and  supersonic  speeds. 

The  effects  of  large  changes  in  principal  axis  inclination  and  mass  dis¬ 
tribution  are  also  included. 

The  primary  purpose  of  this  paper  is  to  summarize  the  information 
obtained  from  the  analog  calculations  and  to  compare  the  results  with 
the  trends  predicted  from  a  slightly  modified  version  of  reference  2. 

SYMBOIE 


b  wing  span,  ft 
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Cz 

Cm 

Cn 

CY 

5 

g 

hP 

it 

XX 

XXe 

XXZ 

XY 

XZ 

M 

M 

m 

N 

P 

P 


lift  coefficient, 

|pv2s 


rolling -moment  coefficient,  Rollins  moment 

|pV2Sb 

pitching-moment  coefficient,  Pitching  moment 

■|PY2S2 


yawing-moment  coefficient. 


lateral-force  coefficient, 


Yawing  moment 
|pV2Sb 

Lateral  force 


|pV2S 


wing  mean  aerodynamic  chord,  ft 
acceleration  due  to  gravity,  ft/sec^ 
pressure  altitude,  ft 
stabilizer  deflection,  deg 

moment  of  Inertia  of  airplane  about  X-axis,  slug-ft^ 

moment  of  inertia  of  rotating  engine  parts  about  X-axis, 
slug-ft^ 

product  of  inertia  referred  to  X-  and  Z-axes,  slug-ft2 
moment  of  inertia  of  airplane  about  Y-axis,  slug-ft2 

moment  of  inertia  of  airplane  about  Z-axis,  slug-ft^ 

Mach  number 

pitching  moment,  ft-lb 

airplane  mass,  slugs 

s 

yawing  moment }  ft-lb 
rolling  velocity,  radians/sec 
average  rolling  velocity,  radians/sec 
steady  rolling  velocity,  radians/sec 
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4 

q 

r 

5 
t 
V 

w 


X,  Y,  Z 


% 

aO 


ctp 

P 

tei,  A0 


e 


4 


p 


pitching  velocity,  radlans/sec 
yawing  velocity,  radians /sec 
wing  area,  sq.  ft 
time,  sec 

true  airspeed,  ft/sec 
airplane  weight,  lb 
body  axes  of  airplane 

angle  of  attack  of  airplane  body  axis,  radians  or  deg 
angle  of  attack  at  zero  lift,  radians  or  deg 

initial  angle  of  attack  of  airplane  body  axis,  deg 

initial  angle  of  attack  of  principal  axis,  deg 

angle  of  sideslip,  radians  or  deg 

increments  from  initial  conditions,  radians  or  deg 

total  aileron  deflection  (positive  for  right  rolls),  radians 
or  deg 

rudder  deflection,  radians  or  deg 

angle  between  body  axis  and  principal  X-axis,  positive  when 
reference  axis  is  above  principal  axis  at  the  nose,  deg 

fraction  of  critical  damping  in  pitch  of  nonrolling  aircraft, 
IfY/.MolY 

fraction  of  critical  damping  in  yaw  of  nonrolling  aircraft, 
-bNr 

mass  density  of  air,  slugs/cu  ft 
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Ap 

% 


<0, 


'e 


°V 


increment  in  angle  of  "bank,  deg 

rotational  velocity  of  engine  rotor,  radians/sec 

nondimens ional  undamped  natural  frequency  in  pitch  of  non¬ 
rolling  aircraft  (ratio  of  pitching  frequency  to  steady 
rolling  velocity) 

nondimens ional  undamped  natural  frequency  in  yaw  of  nonrolling 
aircraft 


Cr  ,  Ct  ,  ,  C-j C7e  j 

W  Lit  H’  *8&t 

Cn^,  Cmp,  Cmg^,  Cinq.^ 
Cnp> 

cYp>  Ma*  NP 

CZp>  clr>  cnp’  cnr> 

CYp>  CYT>  Nr 

Cn>d»  Mq 


V  indicates  derivative  with  respect  to  subscript 


indicates  derivative  with  respect  to 
b 


2V 


x  subscript 


indicates  derivative  with  respect  to 


2V 


X  subscript 


Dot  over  a  symbol  indicates  derivative  with  respect  to  time. 


CALCULATIONS 
Motion  Studies 


The  basic  time  histories,  upon  which  the  analysis  of  this  report  is 
based,  were  calculated  in  real  time  using  a  five-degree-of-freedom  approach 
(forward  speed  assumed  constant).  The  fundamental  equations  used  to  cal¬ 
culate  the  motions  are  shown  in  table  I.  These  equations  were  solved  by 
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means  of  a  Goodyear  Electronic  Differential  Analyzer  ( GEDA) . ^  The  initial 
conditions  used  in  the  calculations  simulated  either  trimmed  level  flight 
or  constant  g  accelerated  flight.  Aileron  deflection  was  generally  the 
sole  control  input  used  in  the  calculations. 

Several  thousand  time  histories  were  run  off  in  the  course  of  the 
general  study.  The  initial  calculations  were  made  with  the  simple  aileron 
input  of  figure  l(a)  (input  A).  A  ramp-type  input  of  50°  per  second  was 
used  until  the  desired  control  angle  was  reached.  At  the  nominal  hank 
angle  (usually  J>60°)  the  control  was  returned  to  zero  at  50°  per  second. 

It  was  found  in  a  number  of  instances  that  this  simple  aileron  input  pro¬ 
duced  unrealistic  results  (for  example,  fig.  1(b))  because  it  failed  to 
arrest  the  rolling  motion  after  the  aileron  was  returned  to  neutral.  To 
obtain  a  more  realistic  evaluation  in  these  instances  it  was  necessary 
to  use  a  control  stick  in  conjunction  with  the  analog,  as  shown  in 
figure  2.  A  typical  control  input  using  the  stick  (input  B)  is  shown 
in  figure  l(c).  It  should  be  noted  that  the  inputs  A  and  B  were  similar 
until  the  point  of  corrective  control.  In  the  latter  instance,  the  oper¬ 
ator  attempted  to  stop  the  rolling  motion  as  rapidly  as  possible. 

A  number  of  approaches  were  tried  in  an  effort  to  arrive  at  the  most 
significant  parameters  for  summary  purposes.  It  was  finally  decided  that 
the  maximum  positive  and  negative  excursions  in  angle  of  attack  and  side¬ 
slip  (fig.  1(c))  plotted  against  the  average  roll  rate  obtained  in  a  given 
maneuver  might  provide  the  best  overall  representation.  The  average  roll 
velocity  is  defined  as  the  bank  angle  at  control  reversal  divided  by  the 
time  required  to  reach  this  angle.  The  bank  angle  was  determined  from 
the  direction  cosine  m^  (table  i). 


Simplified  Analysis 


Throughout  the  analog  program  reference  2  served  as  a  valuable  guide 
notwithstanding  the  fact  that  the  analysis  of  reference  2  assumed  constant 
roll  rate  and  made  other  simplifying  assumptions.  It  was  thought  desir¬ 
able,  however,  to  account  for  engine  gyroscopic  effects  and  to  use  the 


actual  value  for  the  ratio 


instead  of  assuming  unity. 


Another  approach  in  considering  the  effects  of  engine  momentum  on 
rolling  maneuvers  is  presented  in  reference  7» 

It  is  shown  in  appendix  A  that  inclusion  of  the  previously  mentioned 
modifications  to  the  theory  of  reference  2  results  in  an  increment  to  the 

^The  GEDA  equipment  was  made  available  through  the  cooperation  of 
the  Air  Force  Flight  Test  Center. 
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parameters  used  to  define  the  stability  chart,  and  also  lowers  the  hori¬ 
zontal  boundary  slightly  (appendix  A  and  fig.  3).  The  basic  meaning  of 
figure  3  is  Identical  to  that  discussed  in  reference  2  in  that  the  region 
between  a  boundary  and  axis  indicates  a  divergent  condition.  The  roll 
velocity  at  which  the  effective  stability  first  becomes  unstable  for  a 
condition  of  zero  damping  is  referred  to  as  the  resonant  frequency  and 
is  given  by  the  following  equations 


^cnp  |  pV2Sb  IXe(% 

P  =  *  I  '  IX  +  2  (iy  -  IX) 


(1) 


applicable  where  the  vertical  boundary  is  first  intersected  as  the  roll 
velocity  is  progressively  increased,  and 


-I 


-C 


\  pv2sS 


t 


2(lz  -  lx) 


(2) 


applicable  where  the  horizontal  boundary  is  first  intersected. 

Note  that  equations  (l)  and  (2)  are  identical  to  equations  10  and 
11  of  appendix  A}  however,  the  notation  p0  of  reference  2  was  retained 
in  the  appendix  whereas  the  average  roll  velocity  is  used  in  the  remainder 
of  this  paper. 


SCOPE 


Two  basic  fighter  airplane  configurations  are  included  in  this  study. 
One  configuration  had  the  mass  and  aerodynamic  characteristics  represen¬ 
tative  of  current  swept-wing  airplane  types,  while  the  other  used  param¬ 
eters  estimated  for  a  tailless  delta-wing  type.  All  the  rolling  maneuvers 
for  the  swept-wing  airplane  included  in  this  paper  were  calculated  for 
a  Mach  number  of  O.'J,  whereas  the  delta  configuration  is  presented  for 
M  =  0.8  and  M  =  1.2.  Unless  otherwise  indicated,  the  basic  aerodynamic 
and  mass  parameters  shown  in  table  II  were  used  in  the  calculations.  The 
aerodynamic  derivatives  were  based  on  wind-tunnel  results,  flight  tests, 
and  theory.  The  derivatives  used  (with  the  exception  of  Cj  )  were  assumed 


l0‘ 

invariant  with  angle  of  attack.  Table  II  presents  the  variation  of  C} 
with  a  used  for  the  several  basic  flight  conditions. 


0 


For  each  configuration  a  series  of  calculations  were  made  with  suc¬ 
cessive  increments  in  aileron  deflection  using  the  basic  derivatives  of 
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table  II.  Additional  series  were  obtained  by  varying  the  stability  deriva¬ 
tives  Cj^  and  Cnp  and  the  damping  derivatives  and  over 

a  rather  wide  range.  The  effects  of  roll  direction  and  duration  and 
initial  angle  of  attack  were  also  considered  for  each  basic  condition. 


To  obtain  a  better  insight  into  the  mechanism  of  roll  coupling, 
additional  information  was  obtained  for  the  swept-wing  configuration  with 
arbitrary  variations  in  yawing  moment  due  to  rolling  Cnp,  dihedral  effect 

principal  axis  inclination,  altitude,  and  mass  distribution.  The 


importance  of  relatively  small  inadvertent  stabilizer  inputs  was  also 
determined. 


PRESENTATION  OF  RESULTS 


A  summary  of  the  figures  presenting  the  results  of  this  investigation 

are: 


Effects  of: 

Cjj^  at  various  levels  of  Cnp 

Roll  direction  . 

Roll  duration  . 

Damping  in  yaw  . 

Damping  in  pitch  . 

Inadvertent  stabilizer  input  . 
Initial  angle  of  attack  and 
principal  axis  inclination  . 

Altitude  . 

Yaw  due  to  roll  . 

Dihedral  effect  . 

Mass  distribution  . 

Analysis  of  coupling  problem  on 
specific  design  . 


Figures 

4  to 

16 

17  and 

18 

19  to 

21 

22  to 

25 

26  to 

33 

34  and 

35 

36  to 

4o 

4i 

42  and 

43 

44 

45  to 

47 

48 

DISCUSSION 
Individual  Effects 


Importance  of  C^, 


Cnp  relationship.-  Figures  4  to  6  present  a 


summary  of  the  effects  of  variations  in  C^  at  three  levels  of  Cn^ 
for  the  swept-wing  configuration.  The  data  of  figures  4  to  6  present 
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the  envelopes  of  the  maximum  angle  of  sideslip  and  angle  of  attack  excur¬ 
sions  from  trim  (fig.  1(c))  as  a  function  of  average  roll  velocity  in 
nominal  3 6o°  left  rolling  maneuvers.  The  results  indicate  that  the  value 
of  has  a  profound  effect  on  the  calculated  amplitudes  at  all  levels 

of  directional  stability.  Moreover,  the  value  of  at  which  the 


smallest  angles  are  obtained  appears  to  increase  as  Cn^  is  increased. 

The  average  roll  velocity  at  which  the  peak  amplitudes  occur  would  appear 
to  vary  in  a  somewhat  erratic  fashion.  To  analyze  the  trends  shown,  the 
simple  theory  of  reference  2  (modified  as  indicated  in  appendix  A)  was 
utilized.  In  figure  7  the  nondimens ional  frequency  parameters  in  pitch 
and  yaw  are  plotted  with  suitable  terms  added  to  account  for  engine  gyro¬ 
scopic  effects.  For  any  given  flight  condition,  variation  in  the  roll 
velocity  will  trace  out  a  curve  which  for  the  values  of  the  aerodynamics 
and  engine  momentum  considered  is  very  close  to  a  straight  line.  For 
conditions  falling  between  the  coordinate  axes  and  the  vertical  neutral 
stability  boundary  (region  A)  the  simple  theory  predicts  primarily  a 
directional  divergence,  whereas  conditions  between  the  coordinate  axis 
and  the  horizontal  neutral  stability  boundary  (region  B)  are  primarily 
indicative  of  a  pitch  divergence.  It  should  be  noted  that,  inasmuch  as 
the  physical  and  inertia  characteristics  were  assumed  constant,  the  slopes 
of  the  lines  shown  in  figure  7  are  for  all  practical  purposes  dependent 
on  the  ratio  of  to  Cnp.  Each  condition  in  figures  4  to  6  is  num¬ 

bered  and  represented  in  figure  7  by  a  line.  The  several  instances  of 
identical  C ratios  are  indicated  by  a  common  reference  number. 

It  can  now  be  seen  that  the  farther  from  the  origin  the  representative 
line  intersects  the  boundaries  of  figure  7 ,  the  worse  the  coupled  motions. 
Thus,  conditions  1  and  5  which  represent  the  most  extreme  ratios  of  natural 
frequency  in  yaw  and  pitch  also  are  characterized  by  the  least  desirable 
motions.  The  average  roll  velocity  at  which  the  lower  undamped  resonant 
frequency  for  the  nonrolling  airplane  occurs  (equation  1  or  2)  is  shown 
ticked  on  the  curves  in  figures  4  to  6.  It  seems  that  this  frequency  in 
general  occurs  fairly  close  to  the  roll  velocity  at  which  maximum  excur¬ 
sions  occur.  For  a  given  ratio  of  to  Cnn  the  higher  the  roll 


-np 


velocity  for  resonance,  the  greater  the  maximum  amplitude  of  the  motions 
for  the  more  extreme  conditions  2  and  4.  Condition  3>  which  has  a  more 
desirable  proportioning  of  stability,  does  not  show  this  trend  to  any 
degree . 


The  effect  of  increasing  Cn^  at  constant  C^  is  also  evident  in 
figures  4  to  6  and  the  results  follow  the  foregoing  analysis. 

Representative  time  histories  for  aileron  deflections  at  which  near 
peak  amplitudes  were  obtained  are  shown  for  each  combination  of 
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and  Cup  in  figures  8  to  10.  In  nearly  all  instances  the  maximum  excur¬ 
sions  occur  during  the  recovery  phase  of  the  maneuver.  Although  the 
recovery  is  perhaps  not  mathematically  definable,  the  results  of  fig¬ 
ures  4  to  6  were  not  particularly  sensitive  to  the  exact  control  manip¬ 
ulations  used  during  recovery,  and  repeat  runs  produced  very  little 
scatter  on  the  summary  plots.  For  conditions  1  and  2  which  produced  an 
intersection  in  the  upper  portion  of  the  diagram  of  figure  7  the  predom¬ 
inant  angle-of-attack  change  during  the  constant  aileron  portion  of  the 
maneuver  was  negative,  whereas  for  conditions  4  and  5  the  predominant 
angle-of-attack  change  was  positive.  The  initial  sideslip  motion  was 
always  negative,  although  the  yawing  velocity  histories  varied  markedly. 

A  series  of  runs  was  made  in  which  individual  terms  were  eliminated 
singly  from  the  equations  of  table  I.  The  results  are  shown  in  figures  11 
and  12  for  the  two  most  extreme  conditions,  1  and  5  (figs.  4  and  6). 

For  condition  1,  which  represents  a  condition  having  a  Cnp  value 
which  is  much  too  low  in  relation  to  C^,  it  is  evident  that  elimination 

of  the  pq  term  in  the  r  equation  effectively  de-couples  the  motion 
(fig.  11(c)).  When  the  pr  term  in  the  4  equation  is  removed,  however, 
little  effect  is  shown.  The  elimination  of  the  pp  term  in  the  4 
equation  shows  an  effect  similar  to  omission  of  the  pq  term. 

For  condition  5,  which  represents  a  condition  having  a  Cn^  value 
much  too  high  relative  to  C^,  removal  of  the  pr  term  has  the  primary 

de-coupling  effect  (fig.  12).  It  is  seen,  moreover,  that  the  removal  of 
the  pq  or  pP  terms  produces  an  -unfavorable  effect  for  condition  5* 

The  importance  of  the  C^,  Cnp  relationship  for  the  delta-wing 

configuration  at  M  =0.8  and  M  =  1.2  is  illustrated  in  figures  13 
and  l4.  Stability  diagrams  for  these  conditions  are  presented  as  fig¬ 
ures  15  and  l6.  From  a  close  inspection  of  figures  13  to  l6  it  is  seen 
that  the  major  points  brought  out  in  the  discussion  of  the  swept-wing 
configuration  are  substantiated  for  the  delta-wing  configuration. 

For  the  M  =  0.'8  condition  (figs.  13  and  15)  the  expected  favorable 
effect  of  reducing  Cn^  from  -O.36  to  -0.l8  (condition  1  to  condition  3) 

is  evident,  although  the  magnitude  of  the  motions  is  still  quite  large 
for  the  latter  case.  When  the  value  of  is  further  reduced  to  -0.09, 

a  considerable  reduction  in  the  peak  amplitudes  is  obtained  probably 
because  of  the  decrease  in  lower  resonant  frequency. 

At  M  =  1.2,  doubling  the  basic  value  of  Cn^  resulted  in  an  appre¬ 
ciable  improvement  because  of  a  more  desirable  proportioning  of  stability 
(conditions  1  and  3#  fig*  l6)  and  because  of  the  increased  stability  level. 
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Detailed  comparisons  of  the  several  configurations  discussed  should 
be  made  only  in  the  light  of  the  different  initial  conditions  present 
(table  II ).  The  effects  of  some  of  these  variables  will  be  treated  in 
later  sections  of  this  paper. 

Effect  of  roll  direction.-  The  effect  of  roll  direction  is  summarized 
in  figure  17  for  the  swept-ving  configuration  at  Cnp  levels  of  0.057 

and  0.1l4  per  radian.  The  data  indicate  that  in  left  rolls  with  the 
lower  Cnp>  the  peak  motions  occur  at  a  somewhat  lower  average  roll  rate 

and  are  30  percent  larger  than  corresponding  right  rolls.  For  Cnp 

of  0.114  a  similar  trend  is  indicated  but  is  less  marked. 

The  simple  theory  of  reference  2  modified  to  include  the  engine 
terms  has  been  applied  to  the  conditions  of  figure  17.  The  stability 
chart  (fig.  18)  shows  that  for  the  lower  Cn^  the  representative  line 

for  the  left  roll  intersects  in  a  slightly  less  desirable  location  than 
the  corresponding  line  for  the  right  roll.  For  the  higher  CUp  the  two 

lines  are  coincident.  The  lower  resonant  frequencies  for  the  undamped 
conditions  are  indicated  by  the  ticks  (fig.  17)  and  it  is  apparent  that 
the  increase  in  resonant  frequency  calculated  for  the  right  roll-s  cor¬ 
relates  fairly  well  with  the  motion  studies. 

It  should  be  mentioned  that  although  the  left  rolls  of  figure  17 
are  seen  to  be  somewhat  more  critical  than  the  corresponding  right  rolls, 
a  number  of  factors  could  cause  the  reverse  to  be  true.  Among  these 
factors  would  be  direction  of  engine  rotation,  and  the  values  of  CnR  , 

°at 

Ctw  anb  perhaps  initial  angle  of  attack.  For  all  conditions  considered 
P 

in  this  paper  however,  the  left  roll  is  the  more  critical  and  all  of  the 
remaining  results  are  presented  for  left  rolls. 

Effect  of  roll  duration.-  A  summary  of  the  effects  of  roll  duration 
is  shown  for  the  swept-wing  configuration  in  figure  19 .  It  is  seen  that 
the  maximum  amplitudes  calculated  for  the  90°  and  l80°  maneuvers  are  much 
lower  than  for  the  3^0°  rolls.  This  is  not  surprising  if  the  motions  are 
considered  to  be  the  manifestation  of  an  effective  reduction  in  stability 
or  actual  instability,  thus  the  time  duration  (change  in  bank  angle)  would 
be  expected  to  be  a  determining  factor  in  the  motion  build  up.  Represent¬ 
ative  time  histories  are  shown  in  figure  20  for  90°,  l8o°,  360°,  and  l44o° 
maneuvers.  The  possible  effects  of  roll  duration  are  clearly  evident 
particularly  for  the  l44o°  maneuver  in  which  angles  of  attack  and  angles 
of  sideslip  of  very  large  magnitude  are  obtained.  In  most  instances, 
however,  the  amplitudes  attained  in  the  3^0°  rolls  were  comparable  with 
those  calculated  for  l44o°  rolls,  particularly  when  the  recovery  phase 
is  included. 

According  to  simple  theory  (appendix  B)  it  can  be  shown  that  the 
rate  of  divergence  for  an  unstable  condition  will  be  directly  proportioned 
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to  the  roll  velocity  and  also  "be  adversely  affected  when  the  intersection 
of  the  determining  line  with  the  stability  boundary  occutb  at  large  dis¬ 
tances  from  the  "neck"  of  the  diagram.  Thus  the  magnitude  of  the  motions 
might  be  expected  to  build  up  quite  rapidly  for  the  delta-wing  configu¬ 
ration  (condition  1,  fig.  16)  at  supersonic  speeds.  Time  histories 
showing  the  effect  of  bank  angle  are  presented  in  figure  21  for  the  delta¬ 
wing  configuration  at  M  =  1.2.  A  360°  roll  is  seen  to  produce  a  violent 
maneuver  with  sideslip  angles  greater  than  20°  and  large  a  excursions. 
The  l80°  maneuver  still  develops  17-5°  of  sideslip  and  +6°  in  Ax  (±2g). 
Even  the  90°  maneuver,  for  this  condition,  is  accompanied  by  12°  of  side¬ 
slip  but  the  a  amplitudes  are  considerably  reduced. 

Effect  of  damping  in  yaw.-  The  influence  of  changes  in  Cnr  is 
summarized  in  figure  22  for  the  swept-wing  configuration  for  a  Cn^ 
level  of  0.057*  A  tenfold  increase  in  Cnr  reduces  the  peak  p  motions 

appreciably  but  has  very  little  effect  on  the  maximum  angle-of-attack 
amplitudes.  A  value  of  Cnr  100  times  the  basic  value  produced  no  sig¬ 
nificant  additional  decrease  in  the  maximum  (3  amplitudes  although  a 
reduction  in  the  a  motion  is  evident.  Moreover,  particularly  for  the 
highest,  damping  ratio,  the  peak  amplitudes  were  materially  increased  at 
the  lower  roll  rates.  A  stability  diagram  of  the  type  used  previously 
is  shown  as  figure  23*  From  the  simple  theory  it  might  be  expected  that 
a  more  significant  improvement  would  exist  for  the  two  conditions  of 
increased  damping  in  yaw.  Note  also  that  the  peak  amplitudes  of  figure  22 
occur  near  the  same  average  roll  rate,  whereas  the  results  of  figure  23 
might  lead  one  to  believe  the  maximum  excursions  with  increased  damping 
would  occur  at  appreciably  higher  roll  rates. 

Recalling  the  results  of  figure  ll(b),  the  elimination  of  the  rp 
term  in  the  q  equation  did  not  alter  the  basic  level  of  the  motion  for 
conditions  intersecting  the  vertical  divergence  boundary.  For  conditions 
intersecting  the  horizontal  divergence  boundary  (fig.  12(b))  elimination 
of  the  rp  term  was  very  effective. 

Figures  2k  and  25  show  the  effects  of  increasing  Cn^.  by  factors 

of  10  and  100  for  the  extreme  conditions  1  and  5  (figs.  4  and  6).  The 
results  indicate  that  the  increases  in  Cnr  have  a  more  significant 

effect  in  reducing  the  motions  for  condition  5  than  for  condition  1. 

Effect  of  damping  in  pitch.-  The  effects  of  increasing  the  basic 
damping  in  pitch  for  the  swept-wing  configuration  axe  summarized  in  fig¬ 
ure  26  for  two  levels  of  Cn^*  It  is  seen  that  a  threefold  increase  in 
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pitch  damping^  ^9  =  0.27^  results  in  an  appreciable  reduction  in  the 

P  motions  and  some  reduction  in  the  a  amplitudes.  A  tenfold  increase 
in  Cjnq  =  O.92)  produces  a  marked  improvement  at  each  stability- 

level.  The  stability  charts  of  figure  27  show  the  effect  on  the  diver¬ 
gence  boundaries  of  increasing  the  pitch  damping.  These  results  correlate 
fairly  well  with  the  trends  presented  in  figure  26.  The  results  of  the 
time  histories  of  figures  28  and  29  show  the  favorable  effects  of  increases 
in  Cmq  of  10  and  100  times  the  basic  value  on  the  motions  for  the  two 

extreme  basic  conditions  examined  in  the  preceding  paragraph. 

In  figure  30  the  effects  of  variations  in  pitch  damping  are  pre¬ 
sented  for  the  delta-wing  configuration  at  a  Mach  number  of  0.8.  At 
this  Mach  number  increasing  the  level  of  Cm^  showed  a  considerable 

improvement  but  was  not  quite  as  effective  as  for  the  swept-wing  con¬ 
figuration.  The  stability  diagram  is  presented  in  figure  31  for  this 
flight  condition. 

At  a  Mach  number  of  1.2  the  results  shown  in  figure  32  indicate  a 
relatively  small  effect  of  increasing  Cm^  until  an  extremely  large 

level  is  assumed.  It  is  believed  this  condition  is  caused  by  the  inter¬ 
section  of  the  basic  determinant  line  with  the  divergence  boundary  at  a 
relatively  large  value  of  the  ordinate.  It  is  seen  (fig.  33)  that  for 
this  condition  moderate  increases  in  Cm^  would  be  expected  to  have 

only  a  slight  favorable  effect. 

Effect  of  stabilizer  input.-  In  the  previous  section  the  important 
effects  of  pitch  damping  were  discussed.  Inasmuch  as  these  changes  are 
attributable  to  reductions  in  the  pitching  velocity,  it  might  be  reasoned 
that  even  small  stabilizer  inputs  during  the  rolls  could  likewise  have 
an  Important  effect  on  the  results.  Figure  3^  presents  a  summary  of 
results  for  simultaneous  stabilizer  and  aileron  inputs  for  =  0.1l4. 

The  type  of  stabilizer  input  used  is  illustrated  in  the  representative 
time  histories  of  figure  35*  It  is  evident  that  1°  of  airplane  nose-up 
stabilizer  results  in  a  50-percent  increase  in  the  amplitudes  of  the 
(3  motions  and  also  produces  large  increases  in  the  a  excursions, 
whereas  a  stabilizer  input  of  the  opposite  sense  has  a  large  alleviating 


%!he  damping  ratio  considering  only  is  shown  in  the  figures 

because  this  value  was  used  in  calculating  the  boundaries  shown.  Of 
course,  the  total  angle  of  attack  damping  ratio  is  also  dependent  on 

and  the  lift-curve  slope.  Thus,  for  the  swept  wing  airplane  (figs.  26 
and  27)  the  total  damping  ratio  will  be  about  0.l6  higher  than  £0,  whereas 
for  the  delta  wing  configuration  (figs.  30  and.  33)  the  total  damping  ratio 
will  be  about  0.17  higher  at  M  =  0.8  and  0.07  higher  at  M  =  1.2. 


CONFIDENTIAL 


CONFIDENTIAL 


NACA  RM  H56AO6 


Ill- 


effect  (fig.  34).  The  reason  for  these  trends  is,  of  course,  that  the 
negative  stabilizer  increases  the  pitching  velocity  (fig.  35)  while  the 
positive  input  appreciably  reduces  the  pitching  rate.  Reference  8  pre¬ 
sents  flight  data  showing  similar  effects  of  stabilizer  inputs.  When  the 
basic  pitch  damping  is  increased  about  sixfold,  much  less  critical  motions 
are  experienced  for  the  unfavorable  stabilizer  input  (symbols,  fig.  34). 

Data  not  presented  indicated  that  somewhat  more  positive  stabilizer 
inputs  than  shown  in  figure  34  were  undesirable.  This  was  attributed  to 
the  fact  that  large  negative  pitching  rates  were  produced. 

It  should  be  noted  that,  whereas  negative  stabilizer  input  aggra¬ 
vated  the  motions  in  this  instance,  conditions  in  which  the  basic  pitching 
velocity  were  negative  or  instances  in  which  "complementary"  sideslip 
were  present  might  be  made  more  critical  with  positive  stabilizer  inputs 
(table  I). 


It  should  also  be  pointed  out  that  large  changes  in  pitching  moment 
due  to  sideslip  could  also  have  a  decided  effect  on  the  characteristics. 


Effect  of  initial  angle  of  attack.-  The  effects  of  initial  angle  of 
attack  for  the  swept -wing  configuration  are  presented  in  figure  3 6  for 
a  principal  axis  inclination  of  1°  at  two  levels  of  directional  stability. 
Data  are  presented  at  a  =  5°*  1°;  and  “5°  or  at  angles  of  attack  of  the 
principal  axis  of  +4°  and  0°.  The  results  shown  in  figure  36  indicate 
that  the  peak  excursions  occur  at  about  the  same  roll  rate  for  negative 
and  positive  conditions  of  angle  of  attack  with  the  amplitudes  appearing 
as  nearly  mirror  images  of  each  other.  This  is  even  more  evident  in  the 
representative  time  history  plots  of  figures  37(a)  and  37(t>),  which  illus¬ 
trate  the  expected  diametric  opposite  nature  of  the  motions.  For  this 
flight  condition  the  major  disturbing  function  is  produced  by  the  angle 
of  attack  of  the  principal  axis  through  the  Op  term  and  the  I^  terms 


in  the  equations  of  table 


I  ^np  -  Cn5at  -  0  j . 


Thus,  when  the  principal 


axis  is  lined  up  with  the  relative  wind,  essentially  no  motion  is  produced. 
If,  however,  a  small  yawing  moment  is  introduced  by  assuming  Cnp  =  -0.05, 

then  relatively  large  peak^  amplitudes  are  obtained,  particularly  in 
regard  to  sideslip  motions  (symbols,  fig.  3 6).  This  is  not  too  surprising 
if  it  is  considered  that  near  g  =  1.6  radians/sec  an  effective  insta¬ 
bility  is  present  for  the  Cnp  =  0.057  condition.  Thus  a  relatively 


small  out-of-trim  movement  could  be  expected  to  produce  a  sizable  effect. 


^The  yawing  moment  introduced  would  have  produced  about  =  2° 
for  a  roll  rate  of  1.6  radians  if  a  three -degree -of -freedom  uncoupled 
motion  is  assumed. 
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Similar  trends  on  a  somewhat  reduced  scale  are  indicated  for 
Cnp  =  0.114  (fig.  36). 

A  summary  plot  is  presented  in  figure  38  for  a  principal  axis  incli¬ 
nation  of  5°.  The  overall  effects  are,  as  might  he  expected,  with  the 
greatest  amplitudes  recorded  for  conditions  having  the  largest  initial 
angle  of  the  principal  axis.  In  figure  39  the  results  of  figures  36  and 
38  are  compared  by  plotting  the  maximum  values  shown  for  any  conditions 
as  functions  of  the  initial  angle  of  attack  of  the  principal  axis.  The 
importance  of  the  latter  parameter  in  correlating  the  results  is  at  once 
evident  upon  inspection  of  figure  39.  The  amplitudes  were  in  each  case 
a  minimum  when  Op  »  0°. 

It  should  he  mentioned  that  if  rolling  performance  at  initial 
g  levels  of  1  or  greater  is  considered  to  he  of  primary  importance,  then 
a  moderately  large  nose -down  inclination  of  the  principal  axis  (posi¬ 
tive  e)  would  he  beneficial  because  it  would  reduce  Op. 

It  is  believed  that  if  an  appreciable  value  of  Cji  or  Cng  is 

2  at 

present,  the  angle  of  attack  (fig.  39)  at  which  minimum  motions  occur 
would  be  shifted,  depending  on  the  sign  and  magnitude  of  the  additional 
disturbing  functions. 

The  sensitivity  to  change  in  the  initial  angle  of  attack  is  illus¬ 
trated  for  the  delta-wing  configuration  at  a  Mach  number  of  1.2  (fig.  4o). 
Data  are  presented  for  initial  angles  of  attack  of  2.0°  and  3-9° 

(ttp  =  0°  and  1.9°)  •  On  first  inspection  it  might  seem  surprising  that 

such  a  large  motion  would  be  obtained  for  the  case  of  otp  =  0°  at  high 
roll  rate.  However,  a  small  value  of  Cn=  was  used  (table  II )  which 

°at 

provided  the  necessary  disturbing  function.  When  CnR  was  assumed 
zero  (symbols,  fig.  4-0) ,  essentially  no  a  or  p  motions  were  obtained. 

Effect  of  altitude.-  The  influence  of  a  large  altitude  change  on 
the  motions  developed  in  rolls  is  shown  in  figure  4l.  Data  are  presented 
for  the  swept-wing  airplane  (Cn^  =  0.057)  at  initial  angles  of  attack  of 

2°  and  5°  and  altitudes  of  10,000  feet  and  32,000  feet.  It  should  be 
noted  that  the  Initial  angles  of  attack  of  2°  and  5°  correspond  to  level 
flight  conditions  at  10,000  feet  and  32,000  feet,  respectively.  It  is 
evident  from  the  results  that  the  major  effect  of  decreasing  altitude 
(at  constant  initial  a)  is  to  delay  the  buildup  of  large  motions  to  a 
higher  average  roll  velocity.  From  the  loads  standpoint,  of  course, 
much  more  critical  conditions  could  be  obtained  at  the  lower  altitude 
because  of  the  two-  and  one-half -fold  increase  in  dynamic  pressure.  The 
lower  resonant  frequency  calculated  from  equation  (l)  is  noted  by  ticks 
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for  each  altitude  and  it  is  apparent  that  the  increase  in  average  roll 
velocity  at  which  largest  coupling  effects  occur  is  predicted  rather  well 
from  the  simple  theory.  It  might  he  mentioned  that  the  decrease  in  maximum 
amplitude  at  the  lower  altitude  is  possibly  attributable  to  the  60  percent 
increase  in  damping  ratio  in  pitch  /Cmq  assumed  invariant  with  altitude). 


Effect  of  yawing  moment  due  to  roll.-  For  most  of  the  calculations 
a  value  of  Cnp  =  0  was  assumed  for  the  swept-wing  configuration.  Fig¬ 
ure  42  shows  the  effects  of  large  arbitrary  increases  in  Cnp  for 
Cn.  =  0.057  and  0.1l4.  For  the  lower  value  of  Cn_  an  increase  in  Cn 

P  ^p  AJp 

to  0.10  is  seen  to  reduce  appreciably  the  size  of  the  motions  obtained. 

A  further  increase  to  0.20  produces  rather  violent  motions  at  a  relatively 


low  roll  rate.  For  a  value  of 


Cnp  =  0.114 


an  increase  in  C 


np 


to 


either  value  has  a  deleterious  effect  on  the  maximum  amplitudes. 


It  was  mentioned  in  the  discussion  of  initial  angle-of-attack  effects 
that  changes  in  Cj,  or  Cng  which  might  alleviate  the  motions  for 


Jat 


certain  a,Q  could  aggravate  the  coupling  problem  for  other  initial  angles 

of  attack.  Figure  43  shows  the  effects  of  introducing  a  value  of 

cnp  =  0*10  (cnp  =  0.057)  in  rolling  maneuvers  performed  at  Oq  =  5° 

and  -3°.  For  <Xq  =  5°  the  favorable  effect  of  increasing  Cn.p  is 

clearly  shown  (figs.  43(a)  and  (b))  and  probably  results  from  the  large 

change  in  the  yawing  velocity  buildup.  When  the  roll  is  made  at  Oq  =  -3°/ 

however,  the  positive  Cnp  is  seen  to  aggravate  the  motion  appreciably. 


Effect  of  .-  A  number  of  variations  were  investigated  for 

the  basic  swept-wing  configurations  ^Cn^  =  0.057)*  Figure  44(a)  illus¬ 
trates  the  effect  of  doubling  the  ordinates  of  the  variation  with 

a  (table  II).  No  appreciable  change  in  the  amplitude  of  the  motions 
is  indicated. 


Unswept -wing  configurations  are  generally  characterized  by  much 
smaller  variation  in  Ci„  with  a  than  are  swept-wing  configurations. 

It  was  decided,  therefore,  to  study  the  effects  on  the  roll  coupling 
problem  for  three  constant  levels  of  Figure  44(b)  presents  the 

results  for  C7  values  of  0,  -O.063,  and  -0.126  per  radian.  The  data 
*0 

indicate  that  as  the  level  of  C]  is  increased  the  peak  motions  occur 

P 


at  a  slightly  higher  roll  rate,  but  the  peak  amplitudes  are  not  much 
reduced.  Another  aspect  of  the  problem  results  from  the  combination  of 
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high  dihedral  effect  and  adverse  sideslip  which  reduces  the  average  roll 
velocity  attainable  for  a  given  aileron  deflection  (fig.  44).  Thus,  for 
the  highest  Cj  level  if  maximum  aileron  deflection  is  limited  to  30° 


(although  -l4°  of  sideslip  would  be  obtained),  the  angle-of -attack  excur¬ 
sions  would  be  greatly  reduced  at  the  maximum  obtainable  roll  rate  of 
-I.60  radians  per  second  (fig.  44(b)). 


Effect  of  mass  distribution.-  The  importance  of  the  mass  distribu¬ 
tion  parameter  ^T—  is  treated  in  thiB  section.  In  figure  45  data 
J-Z 

are  presented  illustrating  the  effects  obtained  in  changing  from  a 

configuration  with  its  mass  concentrated  along  the  fuselage  5l.-.*x  =  0,9 

Iy  -  IX  IZ 

to  a  condition  of  more  equitable  distribution  — =■= - —  =  0.2.  In  this 

IZ 

study  the  inertia  in  yaw  I^  was  held  constant  (table  III). 


For  a  CnQ  level  of  0.057  it  would  appear  that  reducing  the  mass 
P  Iy  —  Iy 

concentration  along  the  fuselage,  such  that  — ~= - —  is  reduced  from 

0.90  to  0.50,  results  in  a  slight  decrease  in  the  peak  amplitude  of  the 
motions  and  an  increase  in  the  roll  velocity  at  which  these  peak  values 
occur.  A  further  reduction  in  the  inertia  ratio  to  0.20  appears  to 
improve  the  situation  appreciably.  To  compare  these  trends  with  simple 
theory,  a  stability  plot  has  been  prepared  (fig.  46).  The  lower  resonant 
frequency  is  indicated  by  the  ticks  on  figure  4-5.  As  has  been  demon¬ 
strated  previously  in  this  paper,  the  simple  theory  adequately  predicts 


the  trends.  As  the  mass  distribution  factor 


iz 


is  reduced,  the 


resonant  frequency  is  of  course  raised  and  the  intersection  of  the  sta¬ 
bility  line  is  changed  from  one  indicating  an  appreciable  p  divergence 
to  an  intersection  tangent  to  the  a  divergence  side  of  the  boundary. 

A  series  of  representative  time  histories  showing  the  basic  character 
of  the  motions  for  near  peak  coupling  conditions  is  presented  as  fig¬ 
ure  47.  Note  that  the  large  value  of  Ix  associated  with  the  lowest 


inertia  ratio  reduces  the  rolling  acceleration  to  such  a  degree  that  a 
much  higher  aileron  deflection  is  required  to  obtain  a  given  average 
roll  rate  in  the  36O0  maneuvers.  Otherwise,  the  general  character  of 
the  motions  developed  is  quite  similar. 


For  the  Cap  level  of  0.1l4  per  radian,  the  summary  plot  (fig.  45) 

indicates  very  little  difference  in  the  inertia  range  from  0.5  to  0.9. 

Iy  —  Iy 

Moreover,  for  the  condition  of  most  equal  mass  distribution,  — =•= - —  =  0.2, 

XZ 
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the  results  are  seen  to  be  slightly  worse  than  for  the  other  conditions. 
The  reason  for  this  condition  can  he  explained  with  the  aid  of  the  sta¬ 
bility  chart  of  figure  46.  For  this  stability  level,  as  the  mass  is  pro¬ 
gressively  more  evenly  distributed,  the  characteristic  stability  line 
intersects  further  in  the  a  divergence  side  of  the  boundary. 


General  Design  Considerations 

In  the  foregoing  analysis  various  individual  effects  have  been 
examined  relative  to  the  roll  coupling  problem  at  several  flight  condi¬ 
tions.  A  contemporary  swept-wing  airplane  designed  for  Mach  numbers  near 
1.7  is  now  examined  briefly  to  assess  the  roll  coupling  problem  and  also 
to  determine  how  the  results  of  the  previous  sections  apply  to  the  airplane. 

It  should  be  noted  in  the  airplane  design  under  consideration  (tail  C 
of  ref.  8)  the  vertical  tail  size  was  primarily  dictated  by  static  stability 
considerations  at  design  supersonic  Mach  number.  The  longitudinal  stability 
corresponded  to  a  static  margin  of  about  10  percent  at  subcritical  speeds 
with  the  usual  large  increase  in  stability  occurring  at  transonic  speeds. 

First,  a  roll  stability  diagram  is  constructed  (fig.  48(a)),  to  illus¬ 
trate  the  effect  of  Mach  number  on  the  stability  proportioning.  The  deriv¬ 
atives  used  were  taken  from  unpublished  data  (table  IV).  The  stability 
diagram  indicates  that  at  M  =  0.73  the  proportioning  of  longitudinal  to 
directional  stability  is  near  optimum,  whereas  at  M  =  0.93  and  M  =  1.26 
the  normal  rearward  aerodynamic -center  shift  produces  a  less  desirable 
intersection.  The  flight  envelope  of  the  airplane  is  presented  as  fig¬ 
ure  48(b).  The  solid  lines  represent  lines  of  constant  lower  resonant 
frequency  (calculated  from  data  of  table  IV)  and  the  dashed  line  indicates 
the  maximum  capabilities  of  the  ailerons  in  relation  to  the  average  roll 
velocity  obtainable  in  ^>60°  rolls  from  initial  lg  flight. 

For  this  airplane,  in  almost  the  entire  supersonic  flight  range,  the 
average  roll  rates  available  are  considerably  lower  than  the  roll  velocity 
at  which  the  most  violent  maneuvers  might  be  expected.  At  moderate  sub¬ 
sonic  speeds  the  ailerons  are  capable  of  producing  roll  rates  in  which 
serious  roll  coupling  could  be  obtained;  however,  the  desirable  propor¬ 
tioning  of  longitudinal  and  directional  stabilities  would  be  expected  to 
have  an  alleviating  effect  on  the  amplitudes  of  the  motions. 

The  use  of  figure  48  for  a  quick  assessment  of  the  roll  coupling 
problem  is  obvious.  It  should  be  realized,  however,  that  conditions 
falling  slightly  to  the  right  of  the  available  aileron  power  line  can  also 
produce  large  motions  as  indicated  in  previous  sections  of  this  paper. 

This  condition  is  particularly  important  when  the  increase  in  dynamic 
pressure  is  considered.  Preparation  of  a  figure  similar  to  figure  48  at 
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a  number  of  g  levels  supplemented  by  a  knowledge  of  principal  axis  incli¬ 
nation  would  appear  to  be  a  logical  first  step  in  arriving  at  critical 
conditions  for  detailed  analog  calculations. 

As  indicated  in  reference  8,  a  comprehensive  roll  investigation  of 
the  airplane  (full  aileron  560°  rolls)  has  been  made  particularly  at  the 
flight  conditions  shown  by  the  symbols  in  figure  48.  At  supersonic  speeds, 
no  serious  coupling  problem  has  been  encountered  in  flight.  At  M  =  O.73, 
although  no  violent  maneuvers  were  experienced,  the  motions  obtained  were 
fairly  large  and  were  sensitive  to  inadvertent  longitudinal  control  motions. 
At  M  =  0.93  fairly  large  motions  were  also  obtained  and  calculations  indi¬ 
cated  that,  had  slightly  higher  aileron  deflections  been  available,  a  rather 
serious  maneuver  might  have  been  encountered. 

It  would  appear  that  the  airplane  of  reference  8  is  acceptable  from 
the  overall  roll  coupling  standpoint.  Based  on  the  results  of  figures  2 6 
and  34,  however,  a  substantial  increase  in  pitch  damping  would  materially 
improve  the  situation  at  subsonic  speeds  both  for  rolls  made  with  and 
without  inadvertent  stabilizer  inputs.  At  supersonic  speeds  the  aileron 
power  was  considered  more  than  adequate  by  the  pilot,  although  the  average 
roll  rates  were  slightly  lower  than  at  subsonic  speeds. 

To  summarize,  then,  it  would  appear  that  present  design  practice 
(adequate  static  directional  stability  at  design  high-speed  Mach  number 
and  moderate  longitudinal  stability  at  subcritical  speeds )  would  probably 
furnish  a  good  starting  point  in  obtaining  acceptable  motions  in  rolling 
maneuvers.  It  might  be  desirable,  however,  to  insure  the  presence  of  an 
adequate  margin  between  the  lower  resonant  frequency  and  the  average  roll 
velocity  obtainable  from  the  ailerons  at  supersonic  speeds.  In  this  con¬ 
nection  it  should  be  stated  that  the  ailerons  should  be  designed  so  that 
excessive  roll  rates  are  not  attainable  particularly  at  supersonic  speeds. 
Finally,  if  a  pitch  damper  producing  a  total  damping  ratio  on  the  order  of 
0.7  to  1.0  at  subcritical  speeds  is  used,  the  tendency  for  serious  roll 
coupling  effects  in  this  speed  range  will  be  minimized. 


CONCLUSIONS 


A  fairly  comprehensive  five-degree -of -freedom  analog  study  of  the 
roll  coupling  problem  for  a  generalized  swept-wing  and  for  a  tailless 
delta-wing  airplane  configuration  has  been  completed  and  the  following 
conclusions  are  in  order: 

1.  The  various  aerodynamic  and  inertia  parameters  considered  indi¬ 
cated: 

(a)  The  relationship  of  the  longitudinal  and  directional  stabilities 
is  of  paramount  importance.  An  optimum  condition  exists  when  the  resonant 
frequencies  in  pitch  and  yaw  are  approximately  equal. 
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(b)  Increases  in  pitch  damping  had  a  pronounced  favorable  effect  in 
reducing  the  amplitudes  of  the  motions  encountered  and  were,  in  general, 
considerably  more  effective  than  corresponding  changes  in  yaw  damping. 

(c)  Changes  in  mass  distribution  in  the  practical  range  produced 
only  relatively  small  effects. 


(d)  Alterations  in  the  yawing  moment  due  to  rolling  Cnp  although 

sometimes  producing  favorable  effects  would  not  be  useful  in  alleviating 
the  coupling  problem  throughout  the  angle -of -attack  range. 


(e)  Changes  in  the  rolling  moment  due  to  sideslip  parameter 
could  produce  significant  improvements  under  certain  conditions. 


C 


*3 


2.  A  study  of  the  effects  of  flight  condition  revealed  that: 

(a)  The  amplitude  of  the  motions  developed  depends  to  a  large  extent 
on  the  duration  of  the  maneuver.  Limited  studies  indicated  that  90°  roll 
maneuvers  would  be  considerably  less  critical  than  360°  rolls. 

(b)  Smal I  inadvertent  stabilizer  inputs  can  greatly  affect  the  motions 
that  would  otherwise  be  obtained. 


(c)  The  engine  gyroscopic  terms  can  cause  rolls  in  one  direction  to 
be  somewhat  more  critical  than  corresponding  rolls  in  the  opposite  direc¬ 
tion. 


(d)  The  initial  angle  of  attack  of  the  principal  axis  has  an  impor¬ 
tant  bearing  on  the  amplitude  of  the  motions,  particularly  in  the  absence 
of  other  disturbing  functions. 

(e)  For  a  given  initial  angle  of  attack  a  reduction  in  altitude  will 
delay  critical  conditions  to  a  higher  roll  rate,  but  the  maximum  ampli¬ 
tudes  may  be  affected  only  slightly. 

(f )  The  effects  of  variations  in  Mach  number  cannot  easily  be  gener¬ 
alized  since  there  are  likely  to  be  changes  in  most  of  the  controlling 
parameters . 

3.  Utilizing  simple  concepts  proved  useful  in  assessing  the  quali¬ 
tative  effects  of  many  of  the  aerodynamic  and  inertia  parameters  and 
changes  in  flight,  condition.  The  calculated  lower  resonant  frequency 
generally  corresponded  to  the  average  roll  velocity  at  which  the  more 
serious  motions  could  be  expected. 
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4.  It  Is  reasoned  that  rational  design  procedure  can  avoid  the 
problem  of  serious  roll  coupling  at  supersonic  speeds  and  minimize  the 
problem  at  subsonic  speeds  for  the  configurations  of  the  type  considered. 


High-Speed  Flight  Station, 

National  Advisory  Committee  for  Aeronautics, 
Edwards,  Calif.,  December  20,  1955- 
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APPENDIX  A 


If  the  equations  at  the  top  of  page  9  of  NACA  TN  1627  (ref.  2)  are 
written  to  include  the  engine  inertia  terms  and  actual  value  of 

~  — ,  the  equations  become  in  the  notation  of  reference  2 
XY 


0  -  P0t  -  (Poi  +  Po20)5"  +  2 ^0^0  (®  -  Po*)  + 

(t  +  Pq9)  k  +  “02Po20  =  0 


(1) 


where 


*  +  Po0  +  (p02^  -  P0e)F  +  2£^p0  G  +  po0)  - 

(0  -  Pq^)^  +  ^^o2^  =  0 


(2) 


*x  -  % 
iz 


=  F 


IXe“e 

=  K 

% 


Iv  -  I 


X  _ 


Iy 


=  F' 


IXetoe 

iT " Kl 


8 
V 
P  r 


angle  of  pitch  relative  to  flight-path  direction,  radians 
angle  of  yaw  relative  to  flight-path  direction,  radians 
steady  rolling  velocity,  radians/sec 


Then,  as  in  NACA  TN  1627,  assuming 


V  =  p0t  and  D  =  jfr  =  i-  A 


equations  (l)  and  (2)  become 


D2©  -  Di|f  -  D\jrF'  -  0F'  +  ^  +  2(;0a>0D9  -  2£0a>0^  +  cd020  =  0  (3) 


D2*  +  D0  +  fF  -  D9F  -  ^  +  2^Dt  +  2^0  +  cn^  =  0  (4) 
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The  determinant  of  equations  (3)  and  (4)  is 


D2  -  F'  +  2^QtbgD  +  o)fl2  +  JL  -  D  -  DF'  + 


K-jD 

D  -  DF  +  2^  -  -i- 


PO 


Po 


2&^e 


r  1C 

D^  +  F  +  2^D  +  +  -i 


=  0 


P0  =*▼  T"  '  T  p0 

The  determinant  may  he  expanded  to  give  the  quartic 

aD^  '+  bD^  +  cD2  +  dD  +  e  =  0 

The  divergence  boundary  is  obtained  by  setting  the  coefficient  e  =  0 
e  =  -FF'  +  -  F'03t2  +  cn02F  +  H^qCOq  + 


K£  +  »ai£  +  KiK  _  FMCi  +  KjO^e2  =  Q 


P0  p0 


if  ^0  =  ^  =  0  then  equation  (5)  becomes 

e  =  -FF'  +  o)e2ofy2  -  F'Oty2  +  cue 2F  +  —  +  + 

^  ^  PO  Po 

K-iK  F%  Kna),  2 


i  +±0_  =  0 


Po  Po  Po 

from  which  it  can  be  shown 


(5) 


0)Q2  =  F'  -  =  F'  -( 

9  p0  V  XYP0  > 


(6) 


cu,. 2  =  -F  -  ~  =  -F  - 


K- 


■1  _ 


l 

xZP0  , 


(T) 


IXea)e 


if  COg  +  _ 

9  IyPo 


is  plotted  as  ordinate  and  + 


I)^e 

IzPo 


as  abscissa, 
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it  1b  seen  that  a  stability  plot  can  be  obtained  similar  to  those  in 
NACA  TN  1627 


^■Xe^ 
+  ■=■- 

tyPo 


-F 


It  will  be  of  interest  to  determine  the  roll  velocity  at  which  resonance 
occurs.  The  square  of  the  undamped  natural  frequency  in  yaw  can  be 
written 


Cnp  \  PV2Sb 

5 


(8) 


then,  substituting  equation  (8)  into  equation  (7),  enables  calculation 
of  the  steady  roll  velocity  at  which  instability  first  occurs 


Cnp  |  pV2Sb 

W 


=  0 


p . .  t  /ifeLT .  1  ^  (9) 

0  2(Ix  '  XY)  V  -  *1)  \  IX  -  lY  )  9> 
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Inasmuch  as  the  engine  term  inside  the  radical  is  usually  small  relative 
to  the  directional  stability  term,  equation  (9)  can  be  simplified  to 


*>0 

applicable  when  the 
as  roll  velocity  is 


*0 

applicable  when  the 

These  equations  indicate  that  left  rolls  will  have  a  somewhat  lower 
resonant  frequency  than  right  rolls  for  right-hand  engine  rotation. 

For  the  practical  condition  of  finite  damping  in  pitch  and  yaw,  it 
appears  the  boundaries  can  be  calculated  as  if  engine  terms  were  neglected 
(ref.  2)  and  then  used  with  the  modified  parameters  to  account  for  the 
engine  inertia  effects. 


„  j  Cns  1  PV%b 

*  +  2(iy  -  ix) 


IY  -  IX 


(10) 


directional  divergence  boundary  is  intersected  first 
progressively  increased.  Similarly  it  can  be  shown 


/  -Can,  |  pV2S5  IXe% 

TZ  -  XX  +  2(IZ  -  xx) 


(11) 


pitch  divergence  boundary  is  first  crossed. 
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APPENDIX  B 


The  effects  of  roll  rate  and  the  basic  proportioning  of  longitudinal 
to  directional  stability  on  divergence  rate  can  be  examined  by  treating 
a  simple  condition  of  zero  damping  and  by  neglecting  the  engine  terms  of 
equations  1  and  2  of  appendix  A.  If  these  assumptions  sire  made,  the 
determinant  of  equations  3  and  4  of  appendix  A  can  be  expanded  to  give 


iD^  +  cD2  +  e  =  0 


where 


a  =  1 

c  =  1  -  FF'  +  ofy2  +  o>02 
e  =  -FF'  +  oiq2^2  -  o^2F'  +  (d02F 
This  equation  can  be  solved  as  a  quadratic  in  D2  as  follows: 


^2  _  -(Ufy2  +  cue2  +  1  -  FF')  4 
2  - 


hty2  +  o>e2  +  1  -  FF  2 


-  CDg^2  +  co^2F 1  -  0}q2F  +  FF' 


The  value  of  the  positive  real  root  D  of  the  preceding  equation 
will  determine  the  nondimens ional  time  to  double  amplitude  in  accordance 
with  the  formula 


V  =  0-693 

2  Positive  real  root 

where  is  the  nondimens ional  time  to  double  amplitude. 
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Consider  the  specified  points  in  the  unstable  portion  of  the  following 
stability  chart 


If  -F  =  0.71  and  F'  =  0.95  then  the  positive  real  root  for  the 
points  considered  will  be  as  follows: 


Uty2 

0^2 

Positive  real  root 

t'2 

0.5 

2.0 

0.228 

5.03 

.5 

4.0 

•519 

2.17 

.5 

6.0 

_ 

•  557 

1.94 
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Thus  it  would  be  expected  that  if  the  vertical  displacement  of  the 
three  points  considered  is  caused,  for  example,  by  variation  (the 

rate  of  roll  therefore  being  identical)  the  higher  the  location  of  the 
point  in  the  unstable  region  the  greater  the  rate  of  divergence. 

If  on  the  other  hand  the  stability  level  of  both  and  Cnp 

are  altered  such  that  a  given  point  in  the  unstable  region  is  attained 
at  a  different  roll  rate,  it  would  appear  that  the  rate  of  divergence 
would  be  proportioned  to  roll  rate  inasmuch  as 

,  V 

t  2  =  Pcfe  t2  = 

where  t2  is  the  time  to  double  amplitude  measured  in  seconds. 
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ior  of  a  Swept-Wing  Airplane  Having  Lateral-Longitudinal  Coupling. 
NACA  RM  H55L28a,  1956. 
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TABUS  II 

DERIVATIVES  AND  CONSTANTS  REQUIRED  FOR  BASIC  CALCULATIONS 


Swept -wing  airplane 

Delta-wing  airplane 

M  =  0.7 

M  =  0.8 

M  =  1.2 

Basic  flight  condition 

hp,  ft  . 

32,000 

4-0,000 

4o,ooo 

\  pV2,  Ib/sq.  ft  .  . 

197 

175 

395 

V,  ft/sec  . 

690 

777 

1,165 

aQ,  deg . 

5-0 

8.0 

3.9 

ng),  g . 

1.0 

1-33 

1.0 

v  . 

0 

0 

0 

a^,  deg . 

0 

0.8 

1.0 

Physical  characteristics 

c,  ft . 

11.5 

23.13 

b,  ft . 

36.6 

37-8 

S,  sq  ft  . 

377 

661.5 

m,  slugs  . 

745 

745 

I-jj,  slug-ft^  .  .  . 

10,976 

13 

ro 

0 

0 

Iy,  Slug-ft2  .  .  . 

57,100 

106,000 

Iz,  slug-ft2  .  .  . 

64,975 

ll4,6oo 

1X2,  slug-ft2  .  .  . 

942 

3,54-0 

e,  deg  . 

1.0 

2.0 

IXeCJDe>  •  • 

17,55^ 

17,554 
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TABU  II.-  Concluded 


EEttVATIYIS  AND  CONST AlfTB  REQUIRED  FOR  BASIC  CALCUIATIONS 


Swept -ving  airplane 

Delta-wing  airplane 

M  -  0.7 

K  «  0.8 

M  1.2 

Aerodynamic  derivative*  for  basic  flight  conditions 

CLit,  per  radian  ....  . 

0 

0 

0 

Cj^,  per  radian . . . 

3.06 

2.78 

1.8} 

per  radian . 

-1.0 

* 

* 

C^,  per  radian . 

-3-5 

-1.5 

-0.7 

C*^,  per  radian . 

-0.56 

-0.18 

-0.5* 

C^,  per  -radian . 

-1.5 

-0.5 

-0.3 

per  radian . 

0 

0 

0 

per  radian  . . 

0 

0 

0 

C^,  per  radian  . . 

-0.095 

-0.1*0 

-0.150 

Cap,  per  radian  . 

0.05? 

0.070 

0.051 

Cn£,  P er  radian  ....  . 

0 

0 

0 

per  radian . 

0 

-0.01* 

-o.ou 

Ca^,  per  radian . 

* 

» 

* 

CYp,  per  radian  .......  . 

0 

0 

0 

cypJ  P®r  radian . 

0 

0 

0 

Cyp,  per  radian  ............ 

-0.26 

-0.57 

-0.70 

C  yy  per  radian  . . 

0 

0 

0 

per  radian . 

-0.255 

-0.195 

-0.232 

C*r,  per  radian  * . 

0.0*2 

0.070 

0.080 

C2p^j  per  radian . . 

0.05* 

0.0715 

0.050 

C  per  radian  . . . 

* 

* 

» 

Cj^,  per  radian  . 

0 

0 

0 

Cjp,  per  radian  ....  . 

(See  curve) 

(See  curve) 

(See  curve) 

*Kot  applicable 


a»  rudl'inu 

( t»)  IX  ltH-WllK  .  lilrpTUll''. 
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TABLE  III 


MOMENTS  OF  INERTIA  USED  IN  MASS  DISTRIBUTION  STUB! 


IY  -  IX 

Iz 

Iy,  slug-ft2 

Iy,  SlUg-ft2 

Ig,  slug-ft2 

0.20 

27,000 

40,000 

64,975 

.50 

18,200 

51,000 

64,975 

•71 

10,976 

57,100 

64,975 

•  90 

4,430 

63,000 

64,975 

TABLE  IV 


VALUES  OF  Cn  AND  Cm  USED  IN  COMPUTATIONS 
P  "U 


M 

% 

c% 

0.6 

0.100 

-0.39 

.7 

.100 

-.39 

.8 

.100 

-43 

.9 

.110 

-.80 

1.0 

.n4 

-1.14 

1.2 

.103 

-1.14 

1.4 

.063 

-1.00 

Note:  The  physical  and  mass  characteristics  used  to  calculate 
figure  48  were  the  characteristics  of  the  swept-wing 
airplane  of  table  II. 
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p,  radions/sec 


Figure  4. -  Effect  of  Cj^  variation.  Svept-ving  configuration; 

Cnp  =  0.057)  input  B. 
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A/3,  deg  0 


Aa,  deg  of 


.Condition 


'^4  ^8  42  -T6  i25  -24  ^2B  =3 

p,  rodians/sec 

Figure  6.-  Effect  of  C^  variation.  Svept-ving  configuration 

Qnp  =  0*228j  input  B. 
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(a)  Condition  3;  Cb)  Condition  2;  (c)  Condition  1; 

C%=-0.18.  Cma=-0.36.  Cma=-0.T2. 

Figure  8.-  Representative  time  histories  for  swept -wing  configuration 

Input  B;  Cnp  =  0,057. 
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(a)  Condition  4; 


C%  =  -0.18. 


Cb)  Condition  3; 


C!*U  =  -0.36. 


(c)  Condition  2; 


=  -0.72. 


Figure  9.-  Representative  time  histories  for  swept-wing  configuration 


Input  B;  Cnp  =  0.114. 
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(a)  Condition  5>  (b)  Condition  (c)  Condition  3; 

=  -0.18.  c^  =  -0.36.  c,^  =  -0.72. 

Figure  10.-  Representative  time  histories  for  svept-ving  configuration 

Input  Bj  Cnp  =  0.228. 
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Figure  11.-  Relative  importance  of  various  terms  in  equations  of  motions 
Swept-ving  configuration.  Input  Bj  Cj^  =  -0.72j  Cnp  =  0.057* 
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r  equation. 
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a  equation. 


Figure  12.-  Relative  importance  of  various  terms  in  equations  of  motions 
Swept -wing  configurations.  Input  B;  =  -0.l8j  CRp  =  0.228. 
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A/3,  deg 


Figure  14.-  Effect  of  Cj^  and  Cnp  combinations.  Delta-wing 
configuration  M  =  1.2;  input  A. 
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Figure  15. -  Stability  boundaries  for  delta-wing  configurations.  M  =  0.80 

Left  rolls. 
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Figure  18.-  Stability  chart  illustrating  roll  direction  effects.  Swept-wing  configuration. 
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Figure  20.-  Representative  time  histories  shoving  effects  of  roll  duration.  Input  B; 
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C\f/ 

-  -0.095  0.0338 

- -0.95  0.338 


Figure  22.-  Effect  of  increasing  Cnr*  Swept-ving  configuration. 

Input  B;  Cnp  =  0.057* 
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Figure  25,-  Stability  chart  illustrating  effects  of  increasing  Cnr. 
Svept-wing  configuration.  Left  rolls. 
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Figure 


24.-  Time  histories  showing  effect  of  increasing  C 
wing  configuration.  Input  B;  =  -0.72;  Cnp  =  0. 
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Figure  26.-  Effects  of  increasing  Cm  .  Swept-wing  configuration 


Figure  27.-  Stability  charts  illustrating  effects  of  increasing 

configuration.  Left  rolls. 


(a)  -  -3.5-  (*)  -  -35-0.  (c)  - 

Figure  28.-  Time  histories  showing  effects  of  increasing  Cm^. 
wing  configuration.  Input  B;  =  -0.72;  Cnp  =  0.057* 
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Figure  29.-  Time  histories  showing  effects  cof  increasing  C 
wing  configuration.  Input  Bj  C^s  — 0.18;  Cnp  =  0.2 
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-75  0.650 

-15.0  1.320 


Figure  30* -  Effect  of  increasing  C^.  Delta-wing  configuration. 

M  =  0. 80;  input  A. 
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Figure  51.-  Stability  chart  illustrating  effects  of  increasing 
Delta-wing  configuration.  M  =  0.80.  Left  rolls. 
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Figure  33*-  Stability  chart  illustrating  effects  of  increasing 
Delta -wing  configuration.  M  =  1.20.  Left  rolls. 
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p, rodions/sec 

34.-  Effects  of  simultaneous  aileron  and  stabilizer  inputs. 
Swept  -wing  configuration.  Cn^  =  0.114- ;  input  A. 
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Figure  38.-  Effects  of  varying  initial  angle  of  attack.  Swept -wing  configuration 


Figure  4l.-  Effect  of  altitude.  Swept-wing  configuration.  Input  Bj 
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